TOWARDS THE KPP-PROBLEM AND logt-FRONT 
SHIFT FOR HIGHER-ORDER NONLINEAR PDES I. 
BI-HARMONIC AND OTHER PARABOLIC EQUATIONS 



V.A. GALAKTIONOV 

Abstract. The seminal paper by Kolmogorov, Petrovskii, and Piskunov of 1937 [15] 
on the traveUing wave propagation in the reaction-diffusion equation 

(0.1) ut^Uxx + u{l-u) in MxM+, uoix) = H{-x) = {1 ior x < 0; for a; > 0} 

(here H(-) is the Heaviside function), opened a new era in general theory of nonlinear 
PDEs and various applications. This paper became an encyclopedia of deep mathematical 
techniques and tools for nonlinear parabolic equations, which, in last seventy years, were 
further developed in hundreds of papers and in dozens of monographs. 

The KPP paper established the fundamental fact that, in (10. ip . there occurs a travel- 
ling wave f{x — Xot), with the minimal speed Ag = 2, and, in the moving frame with the 
front shift Xf{t) {u{xf{t),t) = i), there is uniform convergence u{xf{t) + y,t) — >• f{y) as 
t +00, where Xf{t) = 2t{l + o(l)). In 1983, by a probabilistic approach, Bramson [3] 
proved that there exists an unbounded log i-shift of the wave front in the PDE problem 
(pTj) and = 2i - I logi(l + o(l)) as t +oo. 

Our goal is to reveal some aspects of KPP-type problems for higher-order semilinear 
parabolic PDEs, including the bi-harmonic equation and the tri-harmonic one 

ut = -Uxxxx + ui'^ ~ u) and Ut ^ u^xxxxx + u{l - u), 

and other poly-harmonic PDEs up to tenth order. Two main questions to study are: 

(i) existence of travelling waves via any analytical/numerical methods, and 

(ii) a formal derivation of the log i-shifting of moving fronts. 



1. Introduction: the classic KPP-problem, known remarkable results, 

and extended pde models 

1.1. The classic KPP— problem of 1937: convergence to TWs. In the seminal 
Kolmogorov-Petrovskii-Piskunov (KPP) problem of 1937 05] 

(1.1) Ut = Uxx + u{l — u) in M X M+, u{x,0) = uo{x) in M, 
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with the step (Heaviside) initial function 

'1, X < 0; 



1.2) uoix) = H{-x) = 



0, X > 0, 



;i-4) 



the solution was proved to converge to the so-called minimal travelling wave (TW) corre- 
sponding to 

(1.3) the minimal speed of TW propagation : Aq = 2. 

Namely, looking for a TW profile f{y) for an arbitrary speed A > yields 

'u*(x, t) = f{y), y = X — Xt, where / solves the ODE problem 
-Af = f" + /(I -/),?/ G M; /(-oo) = 1, /(+oo) = 0. 
This 2nd-order ODE, on the phase-plane {/, /'}, by setting /' = F(/), reduces its order: 

(1.5) ^^ = -x-m^, 

and it was shown that there exists the minimal speed Aq = 2 and the corresponding 
minimalTW profile f{y)- Using the natural normalization 

(1-6) /(O) = i, 

this minimal TW profile is defined uniquely. In addition, 

(1.7) /'(?/)< in R. 

The characteristic equation for the linearized operator in (11. 4p . A = 2, has a multiple zero: 

(1.8) g" + 2g' + g = and = e^^ =^ (;, + 1)^ = o =^ /xi,2 = -1, 
that yields the following asymptotic behaviour of f{y)'- 

(1.9) /(?/) = Co?/ e~^(l + o(l)) as y — )■ +oo, where Cq > is a constant. 

Concerning the relation between the ODE TW problem f 1.5 p for Aq = 2 and the PDE 
Cauchy one (11. ip . (II. 2p . the novel remarkable analysis in [45p of convergence as t — )■ +oo of 



the solution of the Cauchy problem (II. ip . (II. 2p to the minimal TW (II. 4p was performed 
in the TW moving frame. This was very essential, and not in view of the obvious x- 
translational invariance of the equation (II. ip : see below. Eventually, using PDE methods 
and naturally defining the front location via 

(1.10) u{xf{t),t) = l foralM>0, 

the KPP-authors proved that the TW front asymptotically moves like 

(1.11) Xf{t) =2t - g{t), where g'{t)^{} as t ^ +oo. 



^In particular, the original KPP-proof crucially used a version of Sturmian zero set theorem of 1836, 
which was re-discovered by the authors independently. Moreover, in 1937 (101 years after Sturm!), this 
was the first ever use of this Sturm's First Theorem on non-increase of the number of intersections of 
all the TWs for any A > Aq with the solution u{x,t) of (11.11) . (|1.2p of ID parabolic PDEs since Sturm's 
original result in 1836. See further historic comments on those amazing facts in |22[ p. 23] and around. 



Then the convergence result of |15] takes the form: 

(1.12) u{xf{t) + y,t) f{y) as t — )■ +00 uniformly in y G M. 

These classic results were proved for more general source terms Q{u) than the quadratic 
one in f ll.ip . satisfying some natural restrictions. However, here and later on, we keep 
using this simplest reaction term 

(1.13) Q{u)=u{l-u), 

bearing in mind and viewing KPP results as establishing, in particular, a kind of a struc- 
tural stability of the problem. Concerning the ODE (11. 5p . with Aq = 2, the results of [15] 
imply the actual structural stability of the dynamical system (ll.Sp on the phase-plane, 
at least, in a neighbourhood of the heteroclinic connection (0,0) — )■ (1,0), that generates 
the TW profile f{y)- In particular, this means stability relative any small perturbation 
in of the nonlinearity f ll.lSp . A similar much more delicate structural stability result 
is obtained in for the PDE (11. ip relative the source Q. We plan to trace out similar 
(but indeed weaker) structural stability properties for higher-order ODEs involved, so will 
keep dealing with the reaction given by (ll.lSp . 

1.2. Bramson's logt-front drift. We next point out the next principal question arising 
within the KPP ideology. Namely, it is the question on the actual behaviour of the TW 
shift g{t) in (II. lip for t ^ 1, which was not addressed in |15] and, moreover, it was not 
even mentioned therein whether g{t) remains bounded or gets unbounded as t — )■ -|-cxd 
(possibly, it might not that important for the KPP-authors). Later on, it turned out that 
this is a principle, difficult, and rather general question in such PDE problems. Indeed, 
this is about a "centre subspace (manifold) drift" of general solutions of the KPP PDE 
(II. ip along a one-parameter family of exact ODE TWs {/(y + o), a G M}. 

This open problem was solved in 1983 (i.e., 46 years later), when Bramson [3] (see 
also |35]), by pure probabilistic techniques based on the Feynman-Kac integral formula 
together with sample path estimates for Brownian motion, proved that, within the PDE 
setting (II. ip . (II. 2p . there is an unbounded logt-shift of the moving TW frontH 

(1.15) g{t) = k\ogt{l + o{l)) as t ^ +00, with k = |, 

Thus, (I1.15P implies eventual, as t — +00, infinite retarding of the solution u{x^ t) from 
the corresponding minimal TW (uniquely fixed by (II. 6p ). thought the convergence (I1.12p 
takes place in the TW frame. 

^The exponent /c = | in (|1.15p seems deserve to be compared with the "magic" | in Kolmogorov- 
Obukhov theory of local structure of turbulence (discovered in 1941, i.e., just 4 years after the KPP- 
paper). Namely, the famous dimensional Kolmogorov-Ohukhov power "K-Al" law (1941) [44l|49] for the 
energy spectrum of turbulent fluctuations for wave numbers k from the inertial range, 

(1.14) E(k)^Ce^k-i, where e=^{\Du\^), 

describing the rate of dissipation of kinetic energy for high Reynolds numbers Re ((•) is an invariant 
measure of calculating expected values). 

3 



Remark: on PDE approaches to the KPP— 2 problem. Several important results 
have been already proved in many KPP-like studies during last seventy years by using 
PDE methods. For instance, sharp front propagation results have been obtained for 
various classes of positive initial data uo{x) with different decay rates as x — )• +00; see 
[ini [m EH [55] as a source of further references and results. Moreover, it seems, it was not 
studied whether special types of pattern convergence to TWs exist for initial data uo{x) 
of changing sign. For instance, if Uo{x) has a fixed number / > 1 of zeros (sign changing) 
that persist for all times t G M+. In particular, it is unclear still how (bearing in mind 
Sturm's First Theorem of 1836 on nonincrease of the number of zeros; see [22j for the 
history and various applications) could this affect the rate of convergence to the minimal 
TW and the corresponding logt-shift? 

It seems, in this framework, there is no still a fully developed PDE (non-probabilistic 
using ideas of Brownian motion) approach to establishing such a precise expansion as in 

dm]), dni. 

1.3. The main goal of the present paper: extensions of the KPP— problem to 
higher-order semilinear parabolic PDEs. The main goal of the present paper is to 
show that the KPP-ideology is very wide and can be extended (along very similar lines) to 
a variety of other more complicated higher-order semilinear PDEs with the same source- 
type term. It is worth mentioning now, that, since for such PDEs, no Maximum Principle, 
comparison, Sturm's, and other related properties of order-preserving semigroups apply, 
one cannot expect in principle so complete, detailed, and beautiful results as in the classic 
KPP-Bramson study of 1937-1983. 

Of course, there are already many strong results, to be referred to later on, concerning 
a TW-like analysis in Cahn-Hilliard and related higher-order parabolic equations. How- 
ever, our overall goal here and in forthcoming papers is to initiate a more general study, 
and, using any mathematical means (including various analytic, formal, and, inevitably, 
numerical methods), to show that such a general viewing of the KPP ideas makes sense, 
and that many higher-order semi- and quasilinear PDEs of different types inherit some 
deep (but not all) key features of this classic KPP analysis. 

1.4. On some known front propagation results for bi-harmonic and higher- 
order diffusion. Note that, nowadays, there is a vast enough literature devoted to front 
propagation features for fourth-order semilinear parabolic equations such as the Swift- 
Hohenberg one 

(1.16) ut = -{I + D^)^u + au - u\ 

One of the first and detailed such study was performed in 1990 in the monograph [H]. 
Further and more recent papers can be then traced out using the MathSciNet. From 
applications, such KPP-type problems mean that higher-order diffusion is taking into 
account; see references on various fourth- and 2r?T.th-order semilinear and quasilinear par- 
abolic PDEs in [21 ini [151 HSl ISHl IH] . In particular, as further generalizations, in [2j (see 
also references therein and in the current paper), TW profiles were studied for a class of 
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quasilinear thin film equations 

(1.17) ut + {f{u))x = {b{u)u^)x - {c{u)uxxx)x (c(m) > 0). 

Questions of general instabilities of TWs in such semilinear Cahn-Hilliard and other 
related parabolic PDEs (these results are of a special attention in the present study) were 
also already treated in a number of papers; see [551 E21 SS] and references there in. We 
will quote and use later on some other known results. 

However, it seems that the study that is more and directly oriented to the KPP- 
ideology, including possible types of unbounded front shift, was not performed before. 

1.5. Layout of the paper. Thus, we will discuss some aspects of KPP-type problems 
for higher-order semilinear parabolic partial differential equations (PDEs), with the same 
Heaviside initial data. 

In Section [21 as a natural and simplest KPP-like extension, we consider the first model: 
the semilinear bi-harmonic equation, i.e., a fourth-order semilinear heat or a reaction- 
diffusion equation (SHE-4 or RDE-4) 

(1.18) Ut = —Uxxxx + u{l — u) in M X R_|_, u{x,0) = H{—x) in M. 

The corresponding TW with the speed of propagation A is then governed by the following 
fourth-order ODE: 

(1.19) u^ix,t) = fiy), y = x-Xt =^ _A/' = -/W + /(i _ /), 
with the same singular boundary conditions at infinity: 

(1.20) f{y) — )■ and f{y) — )■ 1 as y — )■ ±oo "maximally" exponentially fast. 

This "maximal" decay of f{y) at infinity somehow includes some kind of the remnants of 
a "minimality" of the possible TW profiles, though, for such higher-order equations, any 
direct specification of such a property is difficult to express rigorously (and literally). 

For any A 7^ 0, the problem fll.lQp . fll.20p is of the elliptic type, but it is not a variational 
one. Therefore, we cannot used advanced methods for higher-order ODEs with potential 
operators associated with homotopy-hodograph and other approaches (see key examples in 
[m |50l [53] ) and/or Lusternik-Schnirel'man and fibering theory (see [271128] and references 
therein). Hence, the ODE fll.l9p . though looking rather simple, and, at least, simpler than 
most of related fourth-order ODEs already studied in detail, represents a serious challenge 
and cannot be tackled directly by known tools of modern nonlinear analysis and operator 
theory. 

The KPP-4 ODE problem (ll.lQp . (11.200 on existence of a heteroclinic connection 
(0, 0, 0, 0) I— )■ (1, 0, 0, 0) in a four- dimensional phase space, is more difficult that its second- 
order KPP-2 counterpart (11.40 . so we cannot easily obtain which "minimal" (if any) speed 
Ao occurs for the Heaviside data H{—x). Moreover, we will show that, in the usual sense, 
such a "minimal speed" to be understood in the usual sense, is in fact nonexistent for 
many higher-order equations. Therefore, we first study TW profiles for rather arbitrary 
A G R, and will denote them by 

/ = /(!/; A). 

5 



Here, for the ODE problem fll.lQp . f ll.20p . it is easy to describe directly the exponential 
oscillatory bundles of orbits as y — ?■ ±00. Such an analysis, even more general and 
complete, is well established; see references in Section O Next, by a shooting technique, 
we can justify existence of such TWs for sufficiently small A > 0. It turns out that, for 
such a KPP-4 problem (as well as for many others), the question of existence of a kind 
of a "minimal" speed becomes irrelevant. Moreover, we do not think that, for such PDEs 
and ODEs without any traces of the Maximum and Comparison Principles, the notions 
of "minimal solution/speed" can be properly and so easily defined. Actually, often, these 
essentially do not make sense, so that this part of the classic KPP analysis seems do not 
admit a direct extension. On the contrary, in several cases, numerically, we have observed 
a "maximal" speed Amax > such that 

(1.21) TW profiles f{y; A) exist for all < A < Amax, and nonexistent for A > A^axj 

and rather sharply estimated Amax = Amax("^) > for m = 2,3,4, i.e., up to the eighth- 
order parabolic equation (see fll.24p below). However, we must admit that such a numer- 
ical identification of the Amax is difficult and is not always reliable, since, for higher-order 
ODEs, the infiuence of artificial boundary conditions at the end points of a sufficiently 
large interval of integration cannot be completely eliminated, especially if the solutions 
are highly oscillatory therein. Sometimes, we also cannot guarantee nonexistence of TW 
profiles f{y) for A much larger than Amax? when the bvp4c solver of the MatLab often 
produces a very slow convergence. 

In Section [3l similarly, we study the semilinear tri-harmonic equation (SHE-6): 

(1.22) Ut = Uxxxxxx + u{l — u) in R x IR+, u{x,Q) = H{—x) in M. 

The corresponding TW with the speed of propagation A is then governed by the following 
sixth-order ODE: 

(1.23) n,(x,t) = /(!/), y = x-Xt =^ -A/ = + /(I -/), with ([120]). 

In Section HJ we present some numerical evidence on existence of various TWs and 
there properties for 2mth-order semilinear heat (parabolic) equations (SHE-2m), such as 

(1.24) ut = (-1)'"+^D^'"m + u{l - u), with the ODEs 

(1.25) u.{x,t) = f{y), y = x-Xt =^ _A/' = (-1)-+V(2H + /(i _ 

(plus fll.20p ). and treat higher-order cases for m = 4 and m = 5, i.e., eight- and tenth- 
order heat equations. A related approach based on construction of a majorizing operator 
for 2mth-order parabolic equations for m > 2 as in fll.24p is discussed in Appendix B. 

Finally, in Section [5l we explain how log t-front shift occurs in those higher-order KPP- 
problems. Section [6] is devoted to the omega-limits of PDEs involved. 

6 



1.6. An extension: towards quasilinear parabolic KPP-problems. The present 
paper deals with higher-order semilinear parabolic equations of reaction-diffusion type, 
while, in the second part of this research [2l], we extend some of the above results to the 
quasilinear KPP-477, problem for 

(1.26) ut = -i\u\''u)^^^^ + u{l-u) in R x M+, 
where n > is a parameter, as well as to some other parabolic equations. 

1.7. An extension: to other PDE types and settings. In [25], we will deal with 
higher-order dispersion and hyperbolic equations such as 

(1.27) Uf = —DI^u + u{l — u) and Uu = —Ux:xxx + u{l — u) , etc. 

As for more "exotic" PDE models, as a formal but quite illustrative examples, we consider, 
in ^25J, higher-order dispersion equations and end up with the following one: 

(1.28) D^u = -Dl\ + u{l- u), 
with eleventh-order ODE for the TW profiles 

(1.29) - aV(9) = -/(II) + /(I - /) in M (plus ([I2nD). 

Overall, using the KPP-setting, we refer to such problems as to KPP-(A;,/), where k 
stands for the order of the differential operator in x and / for the order of the derivative 
in t. Therefore, equation fll.28p represents the KPP-(11,9) problem. 

1.8. Summary of Main General Questions. Thus, for several KPP-(/c, Z) problems, 
with A; > 3 and / > 1, the main questions to study, here and in [25] , are mainly: 

(I) The problem of TW existence: existence of travelling waves via any analyti- 
cal/numerical methods and estimating the A- interval of existence: 

(1.30) A = {A G M : there exists a TW profile f{y; A)}. 

Firstly, we prove the positivity, which is easy for all the parabolic models such as (11.241) : 

(1.31) AG A =^ A>0. 

(II) The problem of a "maximal" speed: whether 

(1.32) A is bounded from above? 
In particular, two main questions arise: 

(1.33) is there a maximal speed Amax = sup A > 0? 

Those questions are "remnants" of the KPP setting. It turns out that, unlike in [35], for 
several types of nonlinear PDEs, we have found that 

(1.34) A =(0, Amax), so that Amax = supA, Amax ^ A, and Amm = OGA^A 

(though analytical/numerical proof of nonexistence for all A > Amax requires further 
study, as some other aspects of these problems). 

(III) The logt-SHlFT problem: using general and natural instabilities of such TWs 
f{y; Ao), with a Aq G A, to show and derive the log t-shifting of the moving front in the 
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problem (11. ip . (11 ■2p . connected with a kind of an "(affine) centre subspace behaviour" for 
the rescaled equation. 

Thus, we beheve that the log t-shifting phenomenon is quite a generic property of many 
semilinear or quasilinear 125] KPP-problems of different types. 

(IV) The omega- limit set uj{H): the last and difficult question, which is not being 
properly understood, is the actual behaviour, as t — )■ +oo, of the solution u{x, t) of 
the PDE KPP-problems with data H{—x), or uj{uo), where uq is a bounded "step-like" 
function. Namely, defining, in a natural sense, the w-limit set oj{H) of the properly shifter 
orbit (11.121) . to discuss whether or not 

(1.35) u{H)C{f{.;X), AG A} and/or u{H) = {f{-, Xo), Aq G A}. 

The latter means that the rescaled (shifted) PDE orbit converges to a single TW profile 
with a special "minimal" speed Aq (though, for such higher-order parabolic flows with 
no Maximum Principle and other order-preserving features, a slow evolution within a 
connected subset of different TWs in the former statement should be also taken into 
account). 

Overall, the latter in (11.351) means that, denoting by Es{Xo) the stable subset (a man- 
ifold of stability: all reasonable bounded "step-like" data Uq{x), for which there exists, 
after proper shifting, uniform convergence to /(■; Aq) of some TW profile /(■; Aq) with a 
Ao G A), we want to know whether 

(1.36) H{-x) G EsiXo)-! 

Numerically, this assumes a full-scale numerical experiments in general PDE settings, 
which is not done her^, but it would be rather desirable to be performed by more pro- 
fessional and, possibly, more applied mathematicians in this PDE area. 

2. The basic higher-order KPP-4 problem 

Thus, consider the KPP-(4,1) (or simply KPP-4, that cannot confuse in the parabolic 
case) problem fll.181) . We begin with its ODE counterpart (11.191) . fll.20p . 

2.1. A is always positive. The first simple result is the same as in the KPP-2. 
Proposition 2.1. 

(2.1) If there exists a solution f{y) ^ of ffLTQj) . fOOj) . then A > 0. 

Proof. Multiplying the ODE (11.191) by /', integrating by parts over M, and using the 
boundary conditions (11.201) yield 

.2 2) /(^')' = - / f^'^f + / - = / f'f'" + [4 - ? ] ~Z 

=imr-^-l = -l<o. □ 

■^The author is not an expert in that at all and is not fond of such a research, which could lead to a 
conflict with his current mathematical interests. 
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2.2. TW profiles do exist: numerical evidence. It is convenient to present next 
numerical results, which directly show the global structure of such TW profiles to be, at 
least partially, justified analytically. 

To this end, we used the bvp4c solver of the MatLab with sufficient accuracy and both 
tolerances at least of the order 

(2.3) Tols = 10"^ - 10"^ 

It is important to note that, as the initial data for further iterations, we always took either 
the Heaviside function 

(2.4) initial data for numerical iteration are often H{—y), 

i.e., as in (11. 2p . or its slightly smoother version for a better convergence. This once more 
had to help us to converge to a proper "minimal" profile (indeed, there are many other 
TW profiles), though, of course, this was not guaranteed a priori. We keep this rule for 
all other KPP-(fc, /) problems of in [SU |25], including the hyperbolic ones (I1.27P and more 
higher-order ones such as ( ]1.28p . where initial velocity, acceleration, etc. were taken zero. 

We show the profiles on sufficiently small y-intervals, such as [—100, 200], though, to get 
right singular boundary conditions (ll.20p . we fixed much larger intervals [—600,600], up 
to [—100, 5000], when necessary, and even bigger ones to guarantee a proper convergence 
and revealing TW profiles of the Cauchy problem. 

Figure [H shows TW profiles for sufficiently small A G [0.2, 0.5]. For smaller A = 0.1, see 
below; the results for the smallest A = 0.01 achieved numerically are shown in Figure |9l 
where /(y) gets very oscillatory for ^ 1. 

In Figure El we show TW profiles for larger speeds A G [0.6, 1]. All of them are clearly 
oscillatory as y — > +oo. Concerning the type of convergence to 1 as y — )■ — oo, it looks 
like it is monotone, but more close research below will show its oscillatory character to 
be justified by asymptotic expansion methods later on. 

Both figures above confirm an important property: 

(2.5) the oscillatory behaviour as iy — t- +oo decreases as A t; 

see further comments below. Figure E] shows the character of those oscillations for A = 0.1 
on intervals [0,300] (a) and [0,600] (b). 

Figure m confirms oscillatory convergence to 1 in the opposite limit, as 2/ —oo. While 
(a) still not that convincing and reveals just a single intersection with the equilibrium 
/ = 1, the (b), in the scale 10~^, shows further five intersections with 1. On better 
scales, we can detect a dozen of intersections at least, so the convergence to 1 is clearly 
oscillatory. 

Figure E] shows TW profiles for A larger than 1. The phenomenon of a "maximal" speed 
(jl.2ip was clearly observed, with the following sharp estimate: 

(2.6) 1.27148 < Amax(2) < 1.27149. 

In Figure [6l we show the profiles f{y) for A's close to Amax in (12. 6p . It turns out that, 
for A ~ A~a^x5 ^he TW profiles f{y) remain oscillatory for ?y ^ 1, as Figure |6] confirms. 
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-X f'=-f^^'+f(1 -f): oscillations of TW profiles forX e [0.2,0.5] 




-10 10 20 30 40 



y 

Figure 1. The TW profiles f{y) of (fTTg]) . (HTZm for A = 0.2, 0.3, 0.4, and 0.5. 




-10 -5 5 10 15 



y 

Figure 2. The TW profiles f{y) of ffrT9D . ([L20]) for A = 0.6, 0.7, 0.8, 0.9, 
and 1. 
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■X f'=-fl'''+f(1-f): oscillatory TW profile forX = 0.1 




f(y) X io"^-^^'=-^*''*+^(''-^)- oscillations of TW profile about forX = 0.1 

llllll. 




(a) oscillations on (0,300) (b) oscillations on (0,600) 

Figure 3. Oscillations of the TW profile /(y) for y > 1; A = 0.1. 



-Xf'=-f''''+f(1-f): behaviour of TW profile about 1 forX = 0.1 




f(y) 



-X f'=-f''''+f(1-f): oscillatory TW profile forX = 0.1 



-30 -25 -20 -15 -10 -5 

y 



(a) scale ~ lO'^ (b) scale - 10^^ 

Figure 4. Illustration of an oscillatory convergence of the TW profile f{y) to 
1 for y < -1; A = 0.1. 



Therefore, nonexistence of a /(y) for A = Amax has nothing to do with a possible reduction 
of the dimension (from 3D to 2D) of the asymptotic bundle as y — )■ +00, when two complex 
roots of the characteristic equation become real; cf. Theorem I2.1( ii). 

In Figure m we show the profile /(y) for A = 1.27148 ~ X^^^xj i^)^ its oscillatory 
behaviour for y ~ 35, (b). 

In correlation with the second half of the statement in f ll.2ip . we observed no TW 
profiles /(y) for A > Amax- It should be noted that, since we always take the Heaviside 
function (12. 4p as the initial data for iterations via the bvp4c-solver, the non-convergence 
of the numerical scheme could reflect nonexistence of TWs satisfying the "minimal" re- 
quirement, with a "maximal" decay at infinity, which is available in the singular boundary 
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Figure 5. The TW profiles f{y) of (fTTT^ . (^M) for A G [1, 1.271]. 



conditions f ll.20p . In other words, for A > Amax, there could be other TW profiles f{y) 
with slower decay at infinity. 

2.3. The exponential bundle as y ^ +oo: linearized analysis. In the present 
semilinear [n = 0) case, such an analysis is reasonably easy, though, nowadays, a related 
stability theory using Evans functions have been developed for more complicated nonlinear 
systems; see [2] , with a full list of references therein, with applications to general thin film 
fiows. 

Consider the linearized about zero ODE in (ll.lQp : 

(2.7) - Xg' = -^(^) + g. 

The corresponding characteristic polynomials and equation take the form 

(2.8) g{y) = e^^y =^ /J+(/i,A) = /i^-A^-l = 0. 

The graph of the function H^{fi, A) for various A is shown in Figure [H One can see that it 
keeps a similar form for all A > 0. It is seen that there exists a single real fi = /ii(A) > 
(the unstable mode), a single /i2(A) < (a stable mode), and, for small A > 0, a complex 
stable mode with /i±(A) having Re fi±{X) < 0. 
Indeed, for A = in (12. Sp . we get 

(2.9) i/+(^,0)=/-l = =^ /i±(0) = ±i, 
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X 1 KPP-4: -X r = -f*'*' + f(1 -f), TW profiles for = 1 .2, 1 .271 




Figure 6. The TW profiles f{y) of 01.19p . (11 .20^ is becoming less oscilla- 
tory for y ^ 1 in the limit (12.61) . 



). f'=-f<'"+f(l-f): TW profile for>.=1 .27148 » X 




10 15 



f(y) xio 



-).f=-f''"+f(1-f): TW profile for>.=1 .27148 » X 




(a) f{y) 



(b) oscillations for y ^ 35 



Figure 7. Properties of the TW profile /(y) for A = 1.27148 « A„ 



together with the stable yUi(O) = — 1 and the unstable yU2(0) = 1. The stable (unstable) 
roots for small A > read 

(2.10) ^^(A) = -l-^ + 0(A2) (^,(A) = l + ^ + 0(A2)). 
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Similarly, the computations of complex roots of fl2.8p yield 



(2.11) ^±(A) = -^[1 + 0(A)] ± i[l - ^ + 0(A3)] + 0(A3), Re^±(A) < 0. 



Together with fl2.10p . this creates the whole 3D stable bundle of oscillatory linearized 
orbits, which we have seen in Figures [1] and [2l 

Figure [9] shows that, according to (12. lip , the TW profiles get more and more oscillatory 
when A approaches 0, i.e., when /i±(A) tend to pure imaginary roots in (12. 9p . The graph 
of f{y) for A = 0.01 in Figure |9] was obtained by integration over the interval [—50, 6000]; 
see also further results in the next subsection. 

The above "local" results can be made global by using the following conclusion, that 
differs the KPP-4 problem from the KPP-2 one, for which multiple characteristic zeros 
take place at the minimal speed 2; see (II. 7p . 

Proposition 2.2. For any A G M, the characteristic polynomial (12. 8p does not admit a 
root of the algebraic multiplicity 2. 

Proof. Assume the contrary: there exist two coinciding roots /i G C, so that 

fi'^ — Xfi — 1 = {fi — fiYifJ-"^ + afi + b), with some constants a, 6 G C. 
Equating the coefficients of fx^, y?, fi, and 1 yields the system 

f-2/x + a = 0, 



Obtaining a = 2fi from the first equation and b = from the last one and substituting 
into the second one, we have 



Thus, we claim the following result, which is sufficient for us for further applications: 

Proposition 2.3. (i) At least for all small A > 0, the linearized equation (12. 7p {and hence 
the KPP-4 one in (11.190 for y ^ 1) admits a 3D stable family of oscillatory solutions as 
y — )■ +00, and a ID unstable one of exponentially divergent orbits. 

(ii) For small A < 0, the stable bundle as y ^ +oo is ID, and the unstable one is 3D. 

Using Proposition 12. 2[ one can improve the case (i) to see that it might persist globally 
for m = 2 for all A > 0. However, for some larger m, this might be not the case globally. 
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(2.12) 





|A^ = -2, 



so that then A G M must be complex. □ 



H(h,?l)=h'*-;lH- 1 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



Figure 8. The graphs of the function in (12.81) for various A. 




-10 10 20 30 40 50 



y 

Figure 9. The TW profiles f{y) of ffrTOj) . f OOj) for A = 0.2, 0.1, 0.05, and 0.01. 



KPP-4; - X f'=-f'*'+f(1 -f), TW profile for ).=0.005 





KPP-4: - X f'=-f'*'+f(1 -f), TW profile for 1=0.005 



1 50 200 



(a) on (-20, 250) (b) on (-20, 1500) 

Figure 10. The highly oscillatory TW profile f{y) for A = 0.005. 

2.4. Oscillatory singular limit A — 0^. This has been discussed above. In fact, this 
is also a principle question concerning uniqueness of A-branches of solutions. By careful 
numerical study of the behaviour of f{y; A) as A — t- 0+ we rule out a possibihty to have a 
saddle-node bifurcation at some small As_n > 0, at which, via such a turning point, there 
appear two branches of TW profiles for A > As_n- 

In Figure [TUl we show the TW profile f{y) for A = 0.005, which gets very oscillatory 
for ?/ > 0. We thus claim that: 

there is no a saddle-node bifurcation for small A > and 

(2.13) 

the A-branch ends up as A —i- 0^ in a singular (oscillatory) limit. 

Therefore, in general, nothing contradicts so far the existence of a single A-branch of 
TW profiles. Of course, this well correlates with the nonexistence of a TW at A = 
(Proposition 12. ip and some asymptotic computations above. 



2.5. The 2D stable bundle as y — oo. Then setting / 
yield the following characteristic equation: 



l+g in f ll.igp and linearizing 



(2.14) 



9^'^ + 9 



Therefore, for A 
with the roots 

(2.15) 



and g{y) 
0, we have 2D stable (Re 



=^ //_(//, A) = A/i + 1 = 0. 

> 0) and unstable (Re (■) < 0) bundles 



/i±(0) 



l-i 



±ij and /i±(0) = ±^. 

By continuity, for all small A > 0, there exists 2D stable manifold of the equilibrium 1 
with the roots 



(2.16) 



stable bundle : At+(A) and ;U+(A). 



Unlike Proposition 12. 2[ there exist multiple roots of the polynomial in fl2.14p . that 
follows from fl2.12p with the last equation reading as ... = 1. 
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Proposition 2.4. For the polynomial in fl2.14p . there exist two cases of double roots: 
(2.17) /i± = ±3"^ = ±1.3161... for A± = ±4 ■ 3"^ = ±9.1180... . 



By (12. 6p . the positive value A+ in (I2.17P (at which the bundles may exchange their 
dimensions) plays no role. Both bundles via (I2.15P persists for small A > (A < 0), so: 

Proposition 2.5. At least for all small |A| > 0, the linearized equation (I2.14p {and hence 
the KPP-A one in (11.191) for / ~ 1 for y ^ — 1) admits a 2D stable family of oscillatory 
solutions as y —oo, and a 2D unstable one of exponentially divergent orbits. 

By Proposition 12. 4[ this result is global in A > 0, but this does not help to explain why 
at Amax ("minimal") TW profiles cease to exist. 

2.6. Local blow-up to — oo. In order to verify the global continuation properties of 
stable bundles, one needs to check whether the nonlinear ODE (I1.19P admits blow-up and 
the dimension of such an unstable manifold. To this end, we re-write it down and, as 
usual, neglect the linear lower-order terms, which by standard local interior regularity are 
negligible for |/| ^ 1: 

(2.18) fi^) = -f + f + Xf' = -f{l + o{l)) as /->oo. 
The unperturbed equation has the following exact solution: 

(2.19) f^'^ = -p =^ foiy) = -jyl^ ^ -oo as y^Y,\ 

where ^ is a fixed arbitrary blow-up point. Studying the dimension of its stable 
manifold, we consider a perturbed solution and substitute it into the full equation (12.180 
to get, on some linearization again: 

(2.20) f = fo + e, \e\ « |/| =^ (Fq - y)^^) - 16806 = -840 + ^ + ... . 

One can see that the leading solution consists of balancing the first two and the last term, 
that gives uniquely 

(2-21) e{y) = -'-^^^ + .... 

Hence, in this expansion, no extra arbitrary parameter occurs. Thus, we arrive at: 

Proposition 2.6. For the ODE (11.191) . the stable manifold of blow-up solutions is ID 
depending on a single parameter being their blow-up point 1^ G M. 

2.7. Some general conclusions. The KPP-4 ODE is difficult for an analytic study in 
the corresponding four-dimensional phase-space, so we do not plan here to perform its 
detailed and complete study. 

However, the above local and blow-up analysis allows us to formulate some important 
consequences. Recall that equation (I1.19P comprises analytic nonlinearities only, so by 
classic ODE theory [8], all dependencies on parameters in our problem are given by 
analytic functions. 

17 



(2.22) 



Theorem 2.1. Fix a X E IR+ in the ODE problem f fTTSj) . ffrTOjl .- 

(i) // the dimension of the stable bundle as y ^ +00 is 3D {as in Proposition 12.3( 2)). 
and 2D as y ^ —00 {as in Proposition 12. 5p . then the number of solutions f{y) of the 
ODE problem (11.191) . (11.201) is not more than finite. 

(ii) // the stable dimension of the bundle as y ^ +00 chanqes to 2D, or even becomes 
ID {as in Proposition I2.3( n) for A < 0), then, most probably^, the KPP-A problem does 
not have a solution. 

Proof, (i) Denote by C^t the real coefficients that control the 3D stable bundle at 
y — )■ +00, and by D\^^l the coefficients of the 2D unstable bundle as ?/ — )■ —00. Bearing 
in mind all three Propositions above, to get a global solution f{y), one has to solve the 
following algebraic system: 

'Yo{Cil''^'^) = -00, 

.-Dunst(C'it ' ' ■*) = 0, I = 1,2, 

where the ffist equation means that blow-up does not occur at finite points and hence 
"Fq = —00" is to be understood in a limit sense (the solution remains uniformly bounded). 
Hence, we arrive at a system of three analytic equations with three unknowns C^^' ' . 
Then number of solutions is not more than countable, and the set of solutions is discrete 
with a possible accumulation point at infinity only. 

Moreover, by our blow-up analysis above, it follows that, since the stable modes as 
y — )■ +00 are linearly independent, blow-up always occurs, if these coefficients Cs\' ' 
take sufficiently large values (a stability property of blow-up for quadratic ODE (I2.18p ). 
Therefore, all the solutions of (I2.22p belong to a large "box" in M^, where only a finite 
number of those are available. 

(ii) The result is obvious: then the corresponding analytic system becomes overdeter- 
mined (more equations than the unknowns), so it is inconsistent with the probability 1, 
meaning that its solvability(if any) can be "accidental" only and then is destroyed by 
a. a. arbitrary perturbations of any of the coefficients/functions involved (no structural 
stability). □ 

2.8. Existence: shooting from y = —00. Actually, the system (I2.22p corresponds 
to "shooting from the right to the left-hand side" that was not principal for such an 
algebraic analysis. However, for a complete proof it is easier to shoot from y = —00. 
Indeed, keeping the same notations of "stable" expansion coefficients, we then arrive at 
two algebraic equations only: 

(2.23) 




"^Actually, in a natural sense, with probability zero. Moreover, for A < 0, this means the actual 
nonexistence of TWs, as Proposition 12. II states. 
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where we have used the fact that, aX y = +00, there exists a ID unstable manifold 
governed by the unique coefficient C^nst- Thus, we obtain a system of two equations with 
two unknowns D^^. ' , and such a system admits a rather straightforward analysis on this 
plane. We refer, as a typical example, to [iBi § 4], where a rather similar geometric-type 
analysis applied to a fourth-order ODE arising from a semilinear Cahn-Hilliard equation. 
See also [151 § 7] concerning ODEs in thin film theory. 

The strategy is as follows. For a fixed D^l^ , we use the second parameter D^f^ to satisfy 
the first equation in (12.231) . It is clear that this is always possible, since in both limits 
— )• ±00, the solution f{y; D^^' ) blows up, so that it is global for at least one value 
named D^^\d'^1^). As the next step, we start to vary D^]^ to get rid of the unstable 
component as y ^ +00, i.e., to satisfy the second equation in fl2.23p . In view of a clear 
oscillatory character of all the solutions as y — t- +00, this always can be done. 

Thus, the existence of a proper TW profile only is not a principal difficulty, but unique- 
ness remains a hard open problem still. We do not know any mathematical reason (in 
particular, a kind of a monotonicity property, usually leading to uniqueness), for which it 
must be unique. But numerically, we multiply observed convergence to the same profile 
from variously perturbed initial data as a starting point of nonlinear iterations. 

3. The KPP-(6,1) problem 

We now, more briefiy than above, consider the KPP-(6 ,1) problem (11.221) and its ODE 
counterpart (ll.23p . (I1.20p . Of course, (12. ip remains valid, so we check existence of TWs 
for A > only. 

3.1. Numerical results. Figure [TT] shows TW profiles for A = 0.3, 0.5, and 1. In Figure 
[T2I TW profiles are shown for larger speeds A G [1, 1.35]. All of them are clearly oscillatory 
as ?/ — 7- +00. Next, in Figure [121 we show TW profiles for larger A > 1.4. 

Figure [T^ a) confirms an oscillatory convergence to 1 in the opposite limit, as ?/ — )■ —00. 
Figure (b) explains the same phenomenon as for the KPP-4 (cf. (12. 6p ). but here 

(3.1) m = 3: 2.12110 < Amax(3) < 2.12111... , 

at which TW profiles cease to exist. For A = 2.12110, the TW profile and its oscillations 
for y ^ 1 are shown in Figure [151 

3.2. Dimensional analysis of the bundles and towards existence. This is not es- 
sentially more difficult than in the KPP-4 problem. Namely, consider the linearized about 
zero ODE in ([123]): 

(3.2) -Xg' = g^'^+g. 
The corresponding characteristic polynomial is now 

(3.3) g{y) = e''y =^ //+(^,A)=/i6 + A/x-l = 0. 

The graph of the function H+{fi, A) for various A is shown in Figure [T6l It is quite similar 
to that for the KPP-(4, 1). Moreover, Proposition 12.21 holds for the polynomial in (13. 2p . 
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-?if' = r' + f(1-f): TW profiles for?L = 0.3, 0.5,1 

f(y) 
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Figure 11. The TW profiles f{y) satisfying (OH]), fOOl) for A = 0.3, 0.5, 
and 1. 

Similarly, we study the bundles as ?/ — i- — oo, where there occurs 

(3.4) H+{n,\) = fi^ + Xfi + 1 = 0. 
It admits double roots 

(3.5) /i± = ±5-H = ±0.7647... at A± = ±2 ■ 5"! = ±0.5231... , 

which, again, have nothing to do with the maximal speed in (13. ip . 

Then, for small |A| > 0, the dimensions of stable/unstable manifolds correspond to 
a proper shooting, but with more parameters. Therefore, while we can prove that the 
family of TW profiles is not more than countable (or even finite), the actual proving of 
existence is more difficult for such a multi-dimensional shooting. 

Note also, that the blow-up problem (cf. the ODE (I2.19P ) 

(3.6) = f 

is now more difficult since (13. 6p admits oscillatory solutions, which can be studied as in 
[T5| § 7] by introducing an oscillatory component represented by a periodic function. 

4. Higher-order problems: KPP-(8,1) and KPP-(10,1) 

Those are two examples to illustrate a possibility of such KPP extensions. By (12. ip . 
one needs to consider A > only. 
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X f = f<^' + f(1 -f): TW profiles for?^ = 1 , 1 .1 , 1 .2, 1 .3, 1 .35 
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Figure 12. The TW profiles f{y) of ^HM . fOOl) for A = 1, 1.1, 1.2, 1.3, 
and 1.35. 



= r' + f(1-f): TW profiles aroundX =1 .5692... (bold line) 




xr = r' + f(1-f): oscillatory TW profiles around 1=1 .5692.. 




(a) TW profiles (b) oscillations for y ^ 1 

Figure 13. TW profiles in the KPP-6 problem for A > 1. 

Thus, in the KPP-(8,1), we have the following ODE problem: 

(4.1) (PDE) u, = -Dlu + u{l-u) =^ (ODE) - A/' = + /(I - /). 

Figure [17] shows TW profiles for A = 0.5, 1, and 2, satisfying the ODE (14.11) with the 
singular boundary conditions (ll.20p . The critical existence speed value in the sense of 
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f{y) 



-)Lf = f<'" + f(1-f): f(y)- 1 for ).= 1, 1.1, 1.2, 1.3, 1.35 (scale -10"^) 
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-X f = ('^' + f(1 -f): oscillations of TW profiles about equilibrium 1 



\ I 
\ I 
\ I 



(a) convergence to 1 (b) oscillations for A* = 1.121. 

Figure 14. Oscillations of the TW profile f{y) as y ^ ±oo. 



■X f = r' + f(1 -f): TW profile forX=2.121 10 - X 




_c -If = f'°' + f(1-f): TW profile forX=2.12110 - i 



(^) f{y) (b) oscillations for y ^ 30 

Figure 15. Properties of f{y) for A = 2.12110. 
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([L2ID is now 

(4.2) m = 4 : 2.107 < A„,ax(4) < 2.108. 

Figure ITS] shows f{y) for A = 2.107 from fl4.2p . while Figure [T^ explains its oscillations 
for y ~ 35. 

In the KPP-(10,1), we deal with the following ODE problem: 

(4.3) (PDF) ut = Dl% + uil-u) =^ (ODF) - A/' = /(^^^ + /(I - /). 

Figure |20] shows TW profiles for A = 0.5, 1 and 2, satisfying the ODE (14. 3 p with the 
boundary conditions f ll.20p . 

It is seen that the TW profiles for both problems are close to each other. For instance, 
for A = 0.5, the difference is not more than 2%. This confirms that it is possible to pass 
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KPP-6: about equilibtium 1 , H(n,X)=p.^ + ?l n - 1 
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Figure 16. The graphs of the function in (13.31) for various A = 0, 1, 2, 10. 
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Figure 17. The TW profiles j{y) satisfying gl]), fOOj) for A = 0.5, 1, and 2. 
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Figure 18. The TW profile f{y) satisfying (gl]), (HTZm for A = 2.107. 
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Figure 19. Oscillations of the TW profile f{y) from Figure [18] for A = 2.107. 
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Figure 20. The TW profiles f{y) satisfying (gS]), ffT^ for A = 0.5, 1, and 2. 

to the hmit m — > +oo in the 2mth-order KPP-problem f ll.24p . f ll.25p . and construct an 
analogy of the KPP-problems for an "infinite-order" PDE and ODE. 

A local stability of bundles can be done for both above parabolic models establishing 
for which A's the shooting problem is well posed and the corresponding algebraic systems 
are consistent. Indeed, proving existence, to say nothing of the uniqueness, are open 
problems. 

5. The origin of logt-SHiPT: centre subspace balancing 

In this section, we explain the origin of the log t- shifting (retarding) from the TW for 
the PDE higher-order KPP-problem. 

More correctly, the actual proof of such logt-drift assumes a delicate matching of the 
solution behaviour on compact subsets in the TW variable, i.e., in the Inner Region, 
with a remote Outer Region for y ^ 1, where the influence of the nonlinear term —u^ 
is negligible, and the actual behaviour is governed by the linear bi-harmonic operator. 
However, we do not present here this kind of a procedure of matching of asymptotic 
expansions of those two regions, which is difficult in the present quite general case. 

Thus, we consider a KPP-type problem for a semilinear PDE 



(5.1) 



Ut = Au + u{l — u) in M X ]R+, 
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where A is a proper homogeneous isotropic hnear differential operator satisfying some 
extra conditions specified below. For instance, we may fix, as a simple example, the 
bi-Laplacian operator 

(5.2) A = -Dl (or = D^, or that is the same). 

5.1. Behaviour for small t > 0: first appearance of logt. Given the Heaviside data 
fll.2p . to study fl5.ip . fl5.2p (with m = 2, but any m > 2 fits quite similarly) for small 
t > 0, one needs to use the rescaled independent variables 

(5.3) z = and r = logt— oo as t — t- 0"*". 
Then the rescaled equation takes the form 

(5.4) Vr= {B-\l)+e^v{l-v), where B = -D^ + \zD, + \l . 
Denoting by V{z) the unique stationary solution of the problem 

(5.5) (B-ij)y = in M, y(-oo) = 1, 1/(+oo) = 0, 
and performing the standard linearization in (15. 4 p yield 

(5.6) v = V + w =^ Wr= (B-ll) +e^V{l-V) + 
provided that \w\ <^ 1. Finally, in such an inner region, we obtain 

(5.7) w{z, r) = e^cpiz) + ... =^ (B - i = hiz) ^ -V{z)il - Viz)). 

According to the linearized analysis in Section |2l h{z) has an exponential decay at infinity, 
so, keeping the leading terms, we may assume that 

(5.8) h{z) = _B_e*"^ + ... as z -oo and h{z) = _B+e^+^ + ... as z ^ +oo, 

where B±, b± are some constants and Re6_ > and Re 6+ < 0. 

Since B is known to possess a discrete spectrum in a weighted L^-space [T3j, where 

(5.9) p{z) = 6"°'^'"''^, uq > is a small constant, 

and a complete set of eigenfunctions {i{jk}, the problem for ip in (15. 7p admits the unique 
solution given by the eigenfunction expansion 

(5.10) '/' = I]cfc^fc, where Ck = {K^l), 

{■0^} being the eigenfunction set of generalized Hermite polynomials as eigenfunctions of 
the adjoint operator 

(5.11) W = -Dt-\zD,. 

In f l5.10p . (■,■) is the dual L^-metric, in which (15. lip represents the adjoint to B. By 
(15. Sp . all such integrals converge. Moreover, one can check that 

(5.12) i^{z) = ^e^±^ + ... as z ^ ±oo. 
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Thus, we obtain the following expansion of the solution u{x,t) for the Heaviside initial 
function: 

(5.13) u{x,t) = V{-^r/i) + t'^{jr/i) + ■■■ , in the domain, where \t ip[-^)\ <t: I. 
Hence, the last linearization condition defines the domain of its applicability: by f l5.13p . 

(5.14) \t(^{z)\<t:l =^ |te^±^/*V4| ^ 1 =^ |x| < t3 |logt|. 

Here, we first observe the appearance of the log t- "shift" (now, a factor) in the rescaled 
variable coming from the inner expansion region. However, a further extension of such 
an expansion towards a proper TW-frame (to see how this log t- multiplier is actually 
transformed into the logt-shift) is a difficult procedure, which we cannot treat here. In- 
stead, below, we will catch this log t-shift in a simpler way, starting now directly from the 
TW-frame. 

5.2. TW-frame: linearization and rescaled equation. We now perform a perturba- 
tion analysis in the TW-frame. Then, we need to assume that the corresponding ODE 
problem 

(5.15) _Ao/ = A/ + /(l-/), 

with the conditions (11.201) . admits a unique solution /. Then, attaching the solution 
u{x, t) to the front moving and setting, for convenience, Xf{t) = Xot — g{t), the PDE reads 

(5.16) u{x,t) = v{y,t), y = x — Xot + g(t) =^ Vt = Av + XoVy + v{l — v) — g'(t)vy. 
We next linearize (15.161) by setting 

(5.17) v{y,t) = f{y) + w{y,t), 
that yields the following perturbed equation: 

(5.18) Wt = Btf — g\t)f' — g'{t)wy — nP' , where Bw = Aw + (1 — 2/)!^ + X^Wy. 

Assuming that, in this g{t)-m.ov'm.g frame, there exists the convergence as in fll.l2p . so 
that w{t) — 7- as t — 7- +oo, one can see that the leading non-autonomous perturbation 
in fl5.18p is the second term on the right-hand side. However, as we show, the last two 
terms, though negligible, will define a proper logt-shift of the front. 

5.3. Where logt-shift comes from. First of all, an extra "centre manifold" -like front 
shifting can be connected with the known results on the instability of travelling waves 
in a number of analogous fourth and sixth-order parabolic equations; see [351 ESI |16] 
and references there in. However, of course, this does not directly imply specifically the 
/o^'t-shifting to be justified. In addition, it was always shown that, typically, the corre- 
sponding linearized operators (like B) possess continuous spectra, the fact that essentially 
complicates our further analysis. 

In this connection, note also that the rescaled equation (15.180 is essentially non- autonomous 
in time, so we can hardly use powerful tools of nonlinear semigroup theory; see |47J. How- 
ever, using a formal asymptotic approach, we will trace out some typical features of 
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a centre subspace behaviour, after an extra rescahng and balancing of non-autonomous 
perturbations. 

Thus, as usual (see Introduction), we assume that g'(t) — as t — +00 sufficiently 
fast, i.e., at least algebraically, so that 

(5.19) |/(t)|«|^'(t)| for 

Under the hypothesis (15.191) . the only possible way to balance all the terms therein (in- 
cluding the quadratic one —w"^) for t 1 is to assume the asymptotic separation of 
variables: 

(5.20) w{y,t) = g'{t)ilj{y) + e{t)ip{y) + where \e{t)\ <t: \g'{t)\ as t ^ +00. 
Here, we omit higher-order perturbations. Substituting (15.201) into fl5.20p yields 

g"m + e'{t)^+...= g'mB^-f) 
^ ■ ^ + - i9'it))M^' + - g'mtW + 2y,ij) - e\t)^' + ... . 

Using (I5.19P and (15.201) in balancing first the leading terms of the order 0{g'{t)) yields 
the elliptic equation for ip: 

(5.22) 0{g'{t)) ( = 0(i)) : B^^ - /' = 0. 

Then balancing the rest of the terms in (15.211) requires their asymptotic equivalence, 

(5.23) g"{t) ~ -{g'{t)f ~ e{t), i.e., g{t) = klogt, g'{t) = f , g"{t) = e{t) = ^. 
Then, we obtain the second inhomogeneous singular Sturm-Liouville problem for ip: 

(5.24) O(^) : Bip = ki; + e{ilj' + 

Thus, the first simple asymptotic ODE in (15.231) gives the log t-dependence as in (I1.15p . 
Finally, we arrive at the following system for {ifj,ip}: 

BV' = /', 



(5.25) 



Solving this system, with typical boundary conditions as in (ll.20p . allows then to continue 
the expansion of the solutions of (I5.18P close to an "affine (i.e., shifted via /' on the RHS) 
centre subspace" of B governed by the spectral pair obtained by translation in (I5.15P : 

(5.26) Ao = and Mv) = f'iv)- 
The asymptotic expansion for t ^ 1 then takes the form 

(5.27) w{y,t) = ^ij{y) + ^if{y) + ..., 

which can be easily extended by introducing further terms, with similar inhomogeneous 
Sturm-Liouville problems for the expansion coefficients. 

Final "spectral" remark: on /c-values. Thus, B does not have a discrete spectrum 
(see references at the beginning of Section 15. 3p . so one cannot get a simple algebraic 
equation for k by demanding the standard orthogonality of the right-hand side in the 
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second equation in (15.251) to the adjoint eigenvector i/jq of B* in the L^-metric (in which 
the adjoint operator B* is obtained), hke 

(5.28) k: {ki(j + k'^{tlj' + = 0. 

Therefore, the system fl5.25p cannot itself determine the actual value of k therein. As 
we have mentioned, the latter requires a difficult matching analysis of Inner and Outer 
Regions, which, for the KPP-4 (and all other problems), remains an open problem. 

5.4. On some other related results and references. It is worth mentioning that, in 
similar cases, it is known that, after a suitable time-rescaling and necessary transforma- 
tions, the orbit can approach the center subspace Spanj^/^o} locally in space (in L'^pioc)- 
First results for the semilinear heat equation 

Ut = Au- vF in x M+ {u> 0) 

treated in such a way were proved in SSI ESI IS] , etc. in the middle of the 1980s; a full list 
of references is given in [211 Ch. 2]; see also various chapters in [33]. The first realization of 
the scaling idea in quasilinear degenerate parabolic equations, which is one of the powerful 
tools to study such a stabilization was established by Kamin in 1973, [12]. Sometimes, it is 
called now a renormalization group method. Basic ideas go back to the dimensional ideas 
in nonlinear problems and to the notion of self-similarity of the second kind introduced by 
Ya.B. Zel'dovich at the beginning of the 1950s. See Barenblatt's book [1], where several 
such ideas were discussed first. Of course, it should be noted that rigorous results obtained 
in the above papers are always related to semilinear and quasilinear secon(i-c»r(ier parabolic 
equations and are heavily based on the Maximum Principle. Therefore, any justification 
of such approaches to bi- and poly-harmonic equations is expected to be very difficult. 

One can see that the above elementary conclusion well corresponds to an "(affine) 
centre subspace analysis" of the non-autonomous PDE (I5.18p . and then r = logt naturally 
becomes the corresponding "slow" time variable; see various examples in [3T| [32| 133] of 
such a slow motion along centre subspaces in nonlinear parabolic problems with global 
and blow-up solutions. In the latter case, the slow time variable is 

r = — ln(T — t) — )■ +00 as t approaches finite blow-up time T~ . 

For the semilinear higher-order reaction-absorption equations such as 

(5.29) Ut = —Uxxxx — \u\^~^u in M x R+, with p > 1, 

existence of log t-perturbed global asymptotics was established in [21]. For finite-time 
extinction, with — 1 < p < 1 in (I5.29p . this was done in [23]. For the corresponding 
blow-up problem with the combustion source 

Ut = —Uxxxx + in M X ]R+, with p > 1, 

centre manifold-like log(T — t)-dependent blow-up singularities were constructed in [20] . 
We must admit that a justification of such log t-corrections in any of KPP-2m problems 
with m > 2 is more difficult and has been obtained for global orbits as t — > -|-oo. For 
blow-up orbits, with | log(T — t)\, this remains an open problem. 
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6. When the uj-limit set consists of TW profiles 

Consider the 2mth-order KPP problem for the PDE (11 ■241) . with smooth bounded 
"step-hke" initial data, 

(6.1) u{x,0) = uo{x) in M: tto(x) — )■ 0, x — t- +oo and Mo(a;) — )■ 1, x — )■ — cxd, 
where the convergence is assumed to be sufficiently exponentially fast. 

Theorem 6.1. Let (ll.24p . (16. ip admit a global uniformly bounded solution u{x,t) having 
a "step-like" form for all t ^ 1, i.e., 

(6.2) u{x,t) —> as X — > +00, u{x,t) ^ 1 as x — > — oo 

sufficiently exponentially fast uniformly in t ^ 1. Let, for all t ^ \, the front location 
Xf{t) can be uniquely and smoothly {say, analytically) defined from the equation 

(6.3) u{xf{t),t) = l for t>l, 
with the following asymptotic representation: 

(6.4) Xf{t) = ^ot — g{t), where g'{t) — )■ as t — )■ +cxo, 

where Xq > is a constant. Then Xq G A, i.e., there is a TW profile (maybe, non-unique) 
/(■; Xq) satisfying the ODE (I1.25P . with (ll.20p . and the omega-limit set uj{uq), defined via 
uniform convergence in the TW frame y = x — XqI + g{t), is contained in this connected 
closed family of TW profiles: 

(6.5) u{uo)C{f{.-Xo)}. 

Proof, (i). Given a sequence {tk} +oo, passing to the limit(in a weak sense, implying, 
by parabolic regularity, stronger convergence) in the rescaled perturbed equation (I5.16P 
as {tk + s} — )■ +CXD, one concludes that 

(6.6) v{-,tk + s) ^ uj{-,s) in, say, L;^^(R;R+) (or in Cioc(...), etc.), 

where, by the assumption in (16. 4p . w{y,s) solves (in the classic sense) the corresponding 
unperturbed equation 

(6.7) Ws = Aw + XoWy + w{l - w) in M x M+, 

with data wq G uj{uo). Moreover, by the definition of the front tracking (16. 3p . 

(6.8) w(0,s) = i =^ D^w{0,.s) = for any /c = 1, 2, 3, ... . 

Since the solutions of (16. 6p are analytic in both y and s for s > [18] (see also related 
more recent results in [H]), (16. 8 p implies that w{y,s) is independent of s, so, under the 
above hypothesis, w{y, s) = f{y; Aq). □ 
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Remark: Blow-up in the parabolic problem is possible. The conditions on data 
tto(x) at the beginning of Theorem 16. II are essential. Indeed, if —Uq{x) <^ uq{x) < —1 on 
some interval, then any such solution blows-up in finite time as t — i- G M+: 

(6.9) limt^T- sup^gjj |M(x,t)| = +00. 

Different proofs can be found in [121 129] . Then blow-up can be self-similar [3] and then 
blow-up is not uniform, i.e., 

(6.10) limt^T- infj^eR w(a^, i) = — oo and lim^.^.^- sup^jg^ i/(a;, t) = -|-cxo. 
As an alternative, blow-up may create centre and stable manifold patterns, |20j . 

Remark 1: Is the w-limit set a point? — No answer. For various smooth (or 
analytic) gradient systems, under some extra assumptions, it is known that the omega 
limit sets consists of a single point. We refer to, e.g., as a first such result, to Hale- 
Raugel's approach that applies to guarantee convergence of the orbits of gradient systems; 
see [3^, where a survey and further references are given. This approach essentially relies 
on spectral properties of the linearized operator (main Hale-Raugel's hypothesis is that 
the equilibrium has multiplicity at most one, or > 1 under special hypothesis) and uses 
properties of stable, unstable, and center manifolds, that are difficult to justify for some 
less smooth equations. We also refer to [6], where a similar approach to stabilization is 
used and other references can be found. 

An alternative application is the Lojasiewicz-Simon approach to parabolic equations. 
Namely, the Lojasiewicz-Simon inequality is an effective tool of studying of the stabiliza- 
tion phenomena in various evolution problems. In particular, it completely settles the case 
of analytic nonlinearities; see references in |7]; see also |T71 |39l SOl |5ll |56] and references 
therein. Another approach [30] is based on Zelenyak's ideas (1968) [57] from one di- 
mension that are mainly connected with Lyapunov functions only for gradient dynamical 
systems; see also a short survey and references in [50] . 

Unfortunately, fl6.7p is not a gradient system in the usual and necessary sense (and 
fl5.16p is a not a "perturbed" one, for which there exists a Lyapunov-like functional that 
is "almost" monotone on evolution orbits). In fact, the bi-harmonic equation fll.lSp 
(as well as any of fll.24p ) admits a "pseudo- Lyapunov function", obtained, as usual, by 
multiplication of the equation by ut in L^: 

(6.11) L[u]it) ^ [1 /(i...)2 + / (cix) -lu' + lu')= J{u,r < 0, 

where c{x) is any smooth function with an exponential convergence as x — )■ oo satisfying, 
for the sake of convergence of the second integral in (16.110 . 

(6.12) c(-oo) = I and c{+oo) = 0. 

Unfortunately, L[u]{t) is not bounded below, so cannot be used in the Lyapunov classical 
analysis. This is also easily seen in the TW frame, where, for v = v{y, s), (16. lip reads 

(6.13) ^ [1 J(^^^)2 + J ^c{y + Aot - g{t)) -lu' + lu^)= J{v, - X,Vy - g'{t)vy)^ < 0. 

Observe that, here, c{y + Xot — g(t)) — )■ as t — )■ -|-oo uniformly on compact subsets 
in y. Since such a Lyapunov function has almost nothing to do with TW profiles (i.e., 
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equilibria), we do not think that these inequahties can be of any help in the stabilization 
problem. 

On the other hand, for generalized gradient KPP-problems in Appendix A, where the 
ODEs for / are variational, these approaches can apply; see comments therein. 

Acknowledgement. The author would like to thank D. Williams, as well as other active 
participants of seminars of Reaction-Diffusion and Probability Groups, Department of 
Mathematical Sciences, University of Bath, who stimulated a number of discussions of 
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Appendix A. Some modified variational KPP-problems 

A.l. A more general semilinear bi- harmonic equation. Consider the fohowing bi- 
harmonic equation, with parameters a,/3,7 G M: 

(PDE) : Ut = -U^xxx + OLUxxx + PUxx + lUx + u{l - u), 

^^'^^ (ODE) : A/ ^ -/W + af" + 13 f" + (A + 7)/' = -/(I - /)■ 

Proposition A.l. For a fixed A G M, the ODE in (lA.ip is (formally) variational in the 
metric of the weighted space Ll{R), where 

(A. 2) p = e~"2^ and a G M solves the cubic equation f + /?) = —(A + 7). 

Proof. Using the general form of symmetric (self-adjoint) ordinary differential operators 
we have to find two coefficients pi{y), P2{y), and the weight p{y) such that 

(A.3) Af ^ -1 [{p,f")" + ip2f'y] ^ -i + 2p'J"' + ip'(+P2)f"+p'2f']. 



This yields the following system of linear differential equations: 

El - 
p 

2pi 



r EL — ^ 

p 



(A.4) <^ 



P 



P1+P2 
If = -(A + 7). 

Hence, pi = p, so, the second equation yields the positive weig ht in (jX2l) : 
(A.5) ¥ = =^ p{y) = e-^>0. 

The third equation then defines 
(A.6) p2 = -p'l -I3p = -p" -l3p = -{^ + f3) e-^. 

Finally, the last equation in (]A.4p yields the necessary condition on a: 



(A.7) p'^ = -{X + j)p =^ = ^(A + 7)e-- 
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i.e., comparing with p2 in (jA.6p gives the cubic algebraic equation in ()A.2p . □ 

Thus, expressing all the coefficients in ()A.3P in terms of the weight p given in ()A.2p yields 
the following problem with a symmetric operator: 

(A.8) Af ^ -i [ipf")" + ^ ipf)'] = -/(I - /), 

admitting the (formal) potential: 

(A.9) Hf) = -Upif'r + ^Ipif'? + Ipaf-lf) in H'.nL^. 

On the other hand, for the corresponding rescaled limit (t = +00) parabolic equation in the 
TW- frame (j6.7p . where 

(A.IO) A = A + Ao/ =^ Mf] = m + t IPifT, 

the potential structure of the operator means existence of a good Lyapunov function: 
(All) ^ci>[^](,) = _/p(^,)2<0. 

As was mentioned in Section [6l this allows a further identification of the uj{uq). 
For the former perturbed equation in (j5.16p . similar computations yield 

(A12) = _;p(^^)2 + Ao5'(t) JpVyVt. 

Therefore, a proper control of the second term on the right-hand side is necessary. It perfectly 
converges and is ~ O(p-) as t — > +00 (as the first term), if we trust the formal expansion 
(I5.27p . In this case, <l>[t'](t) is "almost" Lyapunov function that allows to pass to the limit 
and to identify uj{uo) as a point. However, a full rigorous proof of (I5.27P remains an open 
problem. 

A. 2. Application to a particular potential model. We set /3 = 7 = in (|A.ip . so that, 
for a fixed A > 0, we have from ()A.2p 

(A.13) ao(A) = -2A3 < =^ p = e^^""" ^0 as -00. 

The model now reads 

(A.14) ut = -Uxxxx + aQ{\)uxxx + u{l-u) =^ -A/' = -/(^) + ao(A)/'" + /(I - /). 

Then, ()A.13P means that the metric is now suitable for the ODE in ()A.14p with the second 
condition in (|1.20p saying that f{y) — )• 1 as y — > —00. 

Concerning the first condition, /(y) — > as y — > +00, the suitability of the L^-metric (the 
condition / € Lp(]R+)) will depend on the fact whether the necessary asymptotic bundle 
(to be matched with that about the equilibrium 1 as y ^ — 00, as in Section fI7f\ governed 
by the linearized operator, 

(A.15) - ao/' -Xg'-g = Q, g{y) = e'^^ =^ H+{p, A) = / + 2\i^ p^ _ A// - 1 = 0, 

36 



belongs to L^(M+). The characteristic equation in ()A.15p is reduced to a bi-quadratic one: 



(A.16) ^=-iA3+zy =^ 1,4 _ 3 ;^§j,2 _^ _ ;^ ^ Q 

(A.17) =^ ^.^ = 1 Ai±^iAl+l. 

Then > 0, while yields the following "critical" value of the TW speed: 
(A.18) u^=0 =^ A* = = 2.3926.... 

1 

Indeed, at A = A*, the characteristic equation in ()A.15P has a double root = — | j and, 

(A.19) e'^°^ ye'^o^0L2(R+), 

so that, by this functional variational setting, a 2D bundle of exponentially decaying orbits 
as ?/ — > +00 are excluded. Clearly, this diminishes any possibility to have a proper critical 
point of the corresponding functional ()A.9p . as a TW profile. Indeed, the explicit eigenvalues 
(jATTj) allow one to perform further asymptotic analysis, which, together with a similar one 
as y ^ — oo, eventually may describe a precise range of A, for which critical points of $ 
correspond to TWs. But we will not do that and explain our reasons later on. 
Moreover, since the weight ()A.13p is decaying as y ^ — oo, a possible critical point of $(/) 
may have a behaviour therein that has nothing to do with a TW f{y). In other words, in 
this case, a full and a complete classification of all the asymptotics as y ^ — oo of solutions 
of the ODE in ()A.14p is crucially necessary. 

Overall, though potential operators allow to perform a more rigorous analysis of critical 
points, but an extra essential work is inevitable to guarantee the necessary result concerning 
TW profiles. 

However, our numerical experiments show a different and larger than (lA.lSh "critical" value 
of A. Namely, integrating the ODE in ()A.14p on the interval [—100,50] yields a kind of a 
"minimal" speed of existence 

(A. 20) 3.38 < Amin < 3.39, 

such that, for A < Amin, there are no TW profiles f{y). However, the convergence is very 
slow and, actually, a full convergence with Tols~ 10~^ never happened, even using 16000 
points or more. Therefore, we do not know how much the value ()A.20P may depend on 
boundary conditions and the left-hand end point L_ = —100. 

In Figure EH we show TW profiles satisfying (IA.14P for a = 3.39, 4, and 5. Observe a non- 
oscillatory character of solutions in a neighbourhood of the constant equilibria that well 

corresponds to the distribution of the eigenvalues in ()A.17p for A > A,,,. 
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Figure 21. The TW profiles f{y) of (lAlll . (OOll for A = 3.39, 4, and 5. 



A. 3. Non-variational KPP-problem: arbitrary a. Indeed, with a certain extra anal- 
ysis (and certain doubts), the above variational character of the problem ()A.14p . with 
a = ao(A), A > 0, improves its overall mathematical quality. However, we paid so little 
attential to the actual possible rigorous study of (|A.14|) . because the potential functional 
setting here has nothing to do with the existence of TW profiles. 
To show this, consider the same problem but now with an arbitrary a G M: 



(A.21) 



ut 



A/' = -/(^)+ a/'" + /(!-/), 



which is not variational for any a 7^ ao(A) = — 2A3. Note that, for a > 0, the positivity 
result (|2.1|) remains valid and the identity (12. 2 p now reads 



(A.22) 



-A/(/') 



/\2 



-a 



li\2 



< 0. 



It turns out that, numerically, with the full convergence up to 10~^ — 10"''' (unlike the case 
a = ao(A) above), we observe existence of TW profiles for both values a > (Figure [22l 
A = 1) and a < (Figure [23l A = 1). In other words, there exists a local curve of solutions 
originated at a = with existence in some neighbourhood a ~ 0. 

For this A = 1, we have ao{l) = —2, but, for a < 0, we were not able to extend the solution 
for such negative a's. Namely, we found a ''^minimar value 



(A.23) 



0.413 < amin(l) < -0.412 (integrating on [-50, 50]), 

38 



f (y) KPP-4: -X f'= -f''''+a f"'+f(1 -f), TW profiles forX^^ , a=0, 0.4, 0.8 




y 

Figure 22. The TW profiles f{y) of (^KM, (OH]) for A = 1 and a = 0, 0.4, 0.8. 

such that, for a < amin(l)) there is no convergence at all, i.e., / seems to be nonexistent. 
Therefore, here and in \24\ I25j. we always did not follow any "variational temptation" to 
make our analysis more rigorous for some particular and specially designed KPP-problems 
with potential operators, and, instead, concentrated on more general aspects of the KPP- 
ideology and necessary and required related mathematics. 

Appendix B. A linear model: level propagation estimates from above by 

COMPARISON WITH MAJORIZING ORDER-PRESERVING INTEGRAL EVOLUTION 

For the semilinear KPP-2m problem (ll.24p . (jl.2p . there is no the usual comparison of 
solutions for any m > 2. We consider the corresponding linear poly-harmonic equation 

(B.l) ut = i-l)"^+^Dl"'u in R X M+, uo{x) = H{-x), 

which admits certain comparison techniques [29] by constructing a so-called majorizing 
order-preserving operator. This is done as follows. We first present the fundamental solution 
of the poly-harmonic ()B.1|) : 

(B.2) b{x,t)=t-^e'F{y), y = ^. 

Substituting b{x,t) into (jB.ip . one obtains the symmetric (even) profile F as a unique 
solution of the linear ODE 

(B.3) BF = (-ir+iF(2-) + ^F'y+2k^ = M, J^F{y)dy = l. 
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Figure 23. The TW profiles f{y) of (jMH), (OOl) for A = 1 and a = 0, -0.2, -0.4. 

The rescaled kernel F{y) is oscillatory as y ^ oo and satishes a standard pointwise estimate 
|14j . which is convenient to present in the form 

(B.4) \F{y)\ < DoF{y), F{y) = uie-'^olvl" > q in M, where ui = {f^e~'^°^y^" dy)'^ , 

and -Do 5 do are some positive constants depending on m. The normalizing constant wi in 
()B.4p guarantees that the positive kernel F satisfies the standard normalization: 

(B.5) 4F(y)dy = l, 

that will allow us to use it to create below a majorizing order-preserving operator. The 
exponent a is given by 

a = 2^^i € (1, 2) for all m > 2 [a = 2 for m = 1, i.e., for the Gaussian kernel). 

Here Dq > 1 (or Dq — 1 > 0) is called the order deficiency of the poly-harmonic semigroup. 
One can see that D can be made to be equal to 1, D = 1, iS F > 0, i.e., the semigroup is 
order-preserving, and this happens for m = 1 only. 

Next, it is key for us to get a sharp representation of the constant do = do{m). By WKBJ- 
type two-scale asymptotics of the solution F{y) of (1B.3|1 . it is proved that do is associated 
with the root of the algebraic equation 

(B.6) (-l)-a2™-i = 
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It has 2m — 1 different roots {a^}, and in (jB.4p . we choose 
(B.7) —do = maxfc {Rea^ < 0}. 

For the bi-harmonic flow (IB.ip with m = 2, the rescaled kernel estimate (IB.4P is 

(B.8) \F{y)\ < DF{y) = Dujie''^°^y^^^'' in M, where 4 = 3-2-^= 0.2362.... 

In (|B.4p . there appears F as the rescaled kernel of the majorizing order-preserving integral 
evolution equation, 

(B.9) b{x,t) = t^F{y)>0, y = j^, 

whose solutions can be compared with those of (jB.ip . [29j. Obviously, for the linear KPP 
equations (i.e., (jl.24p without the quadratic term —u^) and the corresponding majorizing 
equation, the fundamental solutions are 

(B.IO) h{x,t) = eH{x,t) and h{x,t) = e^{x,t). 

Consider this majorizing equation and the corresponding poly-harmonic flow 

(B.ll) u{t) = b^{t) * uq, uo{x) > and u{t) = * uq for t >0, 

with some initial data no and uq from the corresponding weig 

oo € (0,2(io]) space of correctedness. Subtracting these two and using (IB.4P yield 

u{t) — u{t) = 6* (t) * -Uq - 6* [t) *Uo >K (t) * Uq — \b*{t) \ * \u(j\ 

(B.12) _ _ 

> 6^,(t) *uo - Dob^{t) * \uo\ > e b{t) * (no - Do\uo\), 

whence the following result: 

Proposition B.2. (Comparison) There is comparison of solutions of (jB.lip . i.e., 
(B.13) u{x,t) < u{x,t) in M x for majorizing data uo{x) = Do\uo{x)\ > in M. 

Consider now the majorizing solution 

(B.14) n(x, t) = moe* 6(x, t) corresponding to Dirac's delta function uq{x) = mo6{x), 

where tuq is a sufficiently large constant to be estimated. We now treat ()B.14p as a solution 
of the boundary value problem in x with the symmetry conditions at the boundary 
point X = (of course, at the same time, this is a solution of the Cauchy problem): 

(B.15) u^{0,t)=u^^^{0,t) = ... = Dl"'-\{0,t)=0 for t>0. 

We now choose mo > so large that 
(B.16) Dou{0,t) <u{0,t) for t>0. 

Then the same argument of majorizing comparison guarantees that 

(B.17) u{x,t) <u{x,t) in M+ x M+. 
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This inequality allows to estimate the ^-level propagation p.lOp for u{x,t): indeed, Xf(t) < 
Xf{t), where Xf{t) solves 

(B.18) e4-^F(|^) = i, or F = e'^ ^ as t ^ oo. 

Using ()B.4p . we finally obtain the estimate 

2m~l _ 2m-l 

(B.19) Xf{t)<Xf{t) = d^ t-^d^ logt{l + o{l)), t^oo. 

Therefore, formally (without a proof, since comparison was not performed for full nonlinear 
KPP-problems; see also comments below), we have the following bound for the minimal 
speed: 

(B.20) Ao < = do 2™ and Aq < = S"* • 2" = 2.9512... for m = 2. 

If this minimal speed actually takes plac^, then (IB.19P gives a lower estimate of A: = k{m): 

(B.21) Ao = A,, : k>K = , A; > fc* = 34 • 2^4 = 0.5533... for m = 2. 

Finally, it is worth mentioning that the estimates ()B.20p and ()B.2ip make more sense for 
the KPP-like problems for the semilinear majorizing order-preserving integral, "pseudo- 
differential" evolution equatiorl^ (recall again that it is not a flow, with no a semigroup and 
t-translational invariance available) 

t _ 

(B.22) u{t) = b{t) * no + J h{t - s) * u{s){l - u{s)) ds, 



which is far away from semigroups for the original KPP-problems for higher-order differen- 
tial operators. 

Remark: on comparison in semilinear KPP 2m problem. For m = 1, there is no 
problem, since —v? < 0, so that u{x, t) is always a super-solution. 
Let m > 2, then the integral equation via the Duhamel principle reads 

t 

(B.23) u{t) = e*6(t) * uq - f b{t - s) * {e'u'^{s)) ds. 



Therefore, the final inequality similar to ()B.12p now contains an extra term: 

t 

(B.24) u{t) - u{t) > e* b{t) * (uq - Do\uo\) + J b{t - s) * {e'u'^is)) ds = Ji{t) + J2(t), 



where Ji{t) > by Proposition IB.21 However, since b{x,t) changes sign, in general, we 
cannot control the positive sign of the last term (though > 0) for absolutely arbitrary 
solutions. Thus, any "unconventional" comparison results are not available, and this is well 
known for higher-order parabolic flows. 

we know, by (12. 6p . this cannot happen for m — 2. 
^This kind of equations we meant in Section [Sj thought, it seems that, for the majorizing operator in 
(|B.22[) . the corresponding rescaled operator B in (|5.18l) still does not have a discrete spectrum. 
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But we recall that we need a proper estimate above for t ^ 1, where our solution is assumed 
to be approaching (with a shift) to a TW, i.e., we may assume that 

(B.25) ti(x, t) — 7> /(x — Xf{t)=XQt — g{t), uniformly in x € M as t ^ +00. 

Moreover, we need such an upper estimate on a level set {u{x,t) = ^} for t S> 1. Therefore, 
the second term on the right-hand side of (1B.24P can be arbitrarily sharply approximated 
by 

t 

(B.26) J2{t) ^ fb{t- s)*{e'f{-))ds as t ^ +00, 



and this approximation gets sharper as t increases. Bearing in mind rather smooth TW pro- 
files /(y) (cf. Figures [H [21 O etc.), representing an "almost monotone" (with just slightly 
oscillating tails at y — )• ±00) heteroclinic path of equilibria 1 and 0, and overall "domi- 
nant positivity" of the fundamental solution b{x,t) (J b{x,t)dx = +1), there is a sufficient 
probability that 

(B.27) J2{x,t) > for t ^ 1 and for large enough x, 

at least in an "a. a. sense", in a natural way. If (lR27l) is true, then, by (jR24l) . the front 
propagation estimates ()B.19p - ()B.2ip actually hold. 

The positivity ()B.27p of the integral in (|B.26P deserves further analytic study, though it 
is very difficult, since the TW profiles are not known explicitly. On the other hand, let us 
note that the asymptotic positivity (|B.27p on the front curve {x = Xf(t)} can be studied 
numerically, but, without any doubt, we believe that it actually takes place justifying the 
above TW speed and fc-estimates. 
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